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AN EXPERIMENTAL INVESTIGATION OF TRANSONIC FLOW PAST 
TWO-DIMENSIONAL WEDGE AND CIRCULAR-ARC SECTIONS 
USING A MACH-ZEHNDER INTERFEROMETER ! 


By ArTEUR EARL Bryson, Jr. 


SUMMARY 


Interferometer measurements are given of the flow fields near 
two-dimensional wedge and circular-are sections at zero angle 
of attack at high-aubsonic and low-supersonte velocities. Both 
subsonic flow with local supersonic zone and supersonic flow 
with detached shock wave have been investigated. Pressure 
distributions and drag coefficients as functions of Mach number 
hare been obtained. The wedge data are compared with the 
theoretical work on flow past wedge sections of Guderley and 
Yoshihara, Vincenti and Wagoner, and Cole. 

Jt is shown that the local Mach number at any point on the 
surface of a finite three-dimensional body or an unswept two- 
dimensional body, moving through an infinite fluid, has a 
stationary value at Mach number 1 and, in fact, remains nearly 
constant for a range of speeds below and abore Mach number 1. 
On the basis of this concept and the experimental data, pressure 
distributions and drag coefficients for the wedge and circular-are 
sections are presented throughout the entire transonic range of 


velocities. 
INTRODUCTION 


DIFFICULTIES OF THEORY AND EXPERIMENT IN THE TRANSONIC RANGE 
OF YELOCITIES 

The difficulties inherent in studying transonic flow are 
well-Enown. "Theoretical analysis is made difficult by the 
nonlinearity of the differential] equations of compressible 
fluid motion. This nonlinearity leads to a change-over in 
type of the differential equations from elliptic to hyper- 
bolic when transition is made from subsonic to supersonic 
speeds. Since the essential feature of transonic flow is this 
mixed subsonic-supersonic character, it is obvious that no 
linearization of the differential equations (at least in the 
physical plane) can adequately describe the flow. 

Wind-tunnel studies in the transonic range are made 
difficult bv the large lateral extent of the perturbation flow 
field around bodies in this range. This means that models 
which are small compared with the test section must be 
used. Even then there is still & range of speeds from just 
below 1£,—1 to just above 14, —1 where the model and/or 
its support configuration are “choked,” that is. where local 
supersonic zones embedded in the subsonic field extend from 
the model to the tunnel walls, or, in the supersonic case, 


where embedded subsonic zones extend to the tunnel walls, 
or shock waves, reflected from the walls, impinge on the 
model. Some progress has been made recently in modifying 
wind-tunnel test sections so as to minimize these effects, 
but, on the whole, the majority of good test data in the 
range very close to 14,—1 has so far come from free-flight 
tests. Some good transonic data are available, however, 
from transonic-bump tests made in wind tunnels (reference 
1). Using small models usually results in low Reynolds 
numbers so that difficulty is often experienced in extrapo- 
lating data to full-size Reynolds numbers; this seems to be 
particularly true of the transonic speed range since the ` 
effects of boundary-layer and shock-wave interactions séem 
to be quite large there (references 2 and 3). 

In this paper it is shown that in many instances tests need 
not be made in the region very close to 3f,,=1 since the 
flow in this range can be inferred from testing below and 
&bove this range and using an interpolation based on the 
fact that the local Mach number at any point on the surface 
of unswept two-dimensional bodies and finite three-dimen- 
sional bodies has a stationary value at M, —1. 


EXISTENCE OF POTENTIAL TRANSONIC FLOWS 


Guderley (reference 4) has made a detailed investigation 
of the possibility of smooth transonic flows (i. e., subsonic 
flows with an embedded supersonic zone in which no shock 
waves appear). He proposes that such smooth flows are 
exceptional, that they are discrete cases occurring for only 
particular body shapes at particular free-stream Mach num- 
bers. Any perturbation of the shape with the Mach number 
held constant (or vice versa), Guderley claims, would result 
in a shock appearing in the flow. This bears an analogy to 
the well-known Busemann supersonic biplane which theoret- 
ically has no shocks (and hence no drag) at a discrete value 
of free-stream Mach number and angie of attack (reference 
5, p. 154). Guderley's proposal is still controversial (e. g., 
see the paper of Sears who has made a critical survey of the 
work to date on the existence of transonic potential Hows 
(reference 6)). 


]t is obvious that the potential flow must break down for 
a given body shape at some Mach number less than 1. The 
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argument whether this breakdown occurs precisely when a 
supersonic region first appears on the body or at a slightly 
higher Mach number seems somewhat academic (although 
very interesting), since it is well-known experimentally that 
the drag-rise Mach number (i. e., the Mach number where 
noticeable shocks first appear) is very close to the critical 


Mach number (i. e., the Mach number at which sonic velocity | 


first appears on the body) for most bodies without surface 
slope discontinuities. 

Kuo (reference 7) proposes that supersonic compression 
is unstable to disturbances; that is, & supersonic region on 
a body in subsonic flow must end in a shock with no com- 
pression occurring in the supersonic flow ahead of the shock. 
There seems to be ample experimental evidence to show 
that this is not strictly true since, for example, the com- 
pression region of a A-shock is clearly supersonic. However, 
the A-shock configuration is believed to be a phenomenon 
associated with laminar-boundary-layer and shock-wave in- 


teraction; with turbulent boundary layer (a condition more ` 


closely approaching nonviscous flow) hardly any noticeable 
supersonic compression occurs before the shock ending the 
supersonic zone (see reference 8). 


CHOICE OF MODELS 


Two-dimensional flow is much simpler to handle than 
axially symmetric flow both in theoretical work and in inter- 
ferometry. Hence it was decided to study two-dimensional 
flows despite the well-known difficulties in approximating 
two-dimensional flow in a wind tunnel. 

Because of the considerations mentioned previously it 
was decided to test very small models which would be of 
such a shape that viscous influences would not materially 
affect the flow over them. This led to the choice of “half 
airfoils’—wedges and circular-arc sections followed by 


straight sections. These models have favorable pressure. 


gradients on their surfaces over most of the transonic 
range so that boundary-layer separation, if it does occur, 
will only occur because of shock-wave influence. Further- 
more, such separation will occur downstream of the part of 


the body being studied and hence will not affect the measure- 


ments. Certain viscous effects will still be evident, however, 
for instance, the effective rounding off of the shoulders and 
leading edges of the wedge models. $ 

Both theoretical advantages and practical need make the 
study of thin sections desirable. Consequently, the semi- 
wedge angles chosen were 44°, 74°, and 10° (a 26.6° wedge 
was also used in order to make a comparison with some avail- 
able theoretical work on a wedge of thisangle). The circular- 
are section chosen was essentially the front half of an 8.8- 
percent-thick biconvex circular-arc airfoil, followed by a 
straight section. Models of sections much thinner than this, 
with the same chord lengths used, run into structural 
difficulties and also the ratia of boundary-layer thickness to 
model thickness becomes large enough to cause considerable 
deviation from nonviscous flow. 


TRANSONIC-FLOW THEORY AND EXPERIMENTS 


The investigations of Von Kármán, Busemann, Guderley, 
Frankl, and many others have contributed significantly to 
methods of approach which can be used to study transonic 
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flow (references 9 to 14). The detailed numerical calculations 


for specific cases made by Maccoll and Codd, Emmons, ` 


Drougge, Drebinger, Guderley and Yoshihara, and Vincenti 
and Wagoner (references 15 to 21) have helped to dispel the 
idea of a “sonic barrier.” Recently Cole at GALCIT has 
given an analysis of the flow past wedge sections at high- 
subsonic speeds (reference 22). By combining the results of 
Guderley and Yoshibara’s, Vincenti and Wagoner's, and 
Cole’s calculations, the flow past thin wedge sections can 
be given completely through the transonic range permitting 
& comparison with the present experiments. Some of the 


investigations mentioned above will be discussed in more 


detail further on in the present paper. 

Available experiments in the transonic range on thin wedge 
sections are surprisingly few. Pack (reference 23) describes 
some interferometric experiments on 10° and 20° semiangle 


wedges made at Braunschweig. His subsonic data appear | 


to be good, but the flow in the supersonic interferograms 
appears to be somewhat nonuniform and not very closely 
two-dimensional; only one supersonic Mach number was 
tested where detached shocks occurred. His conclusion that 
the p/p, distributions on the surface of the 20? semiangle 
wedge are very much the same for Mf,=0.803 and AL zs 
1.40 is interesting, but the statement that this agrees with 
the theoretical predictions of Maccoll and Codd is incorrect 
since they indicated that the p/p, distributions would be 
nearly the same. 

Griffith at Princeton has just recently published tho 
results of some very carefully done experiments on flow past 
wedge sections of semiangles of 7°, 10°, 20°, 30°, 45°, and 
90° (and several other shapes) with detached shock waves 
(reference 24). These experiments were,done in a shock 


tube and interferograms are presented of the flow fields. ` 


The experiments clearly show that the shape of the detached 
shock and its detachment distance from the sonic point on 
a wedge depend only on the body thickness and the Mach 
number (not the wedge angle) when the Mach number is 
well below the shock-attachment Mach number. This is in 
general agreement with Busemann’s 
paper on detached shock waves (reference 10). 

Liepmann, Ashkenas, and Cole (reference 8) made some 
careful pressure measurements on the surfaces of 6- and 
12-percent-thiek biconvex circular-are airfoils at zero angle 
of attack at high-subsonic speeds in connection with studies 
of shock-wave and boundary-layer interaetion. Some of the 
results of their tests are combined here with corresponding 
low-supersonic test results from the present investigation to 


indicate the behavior of the pressure distribution on circular- | 


arc airfoils at zero angle of attack through the entire transonie 
range. 
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SYMBOLS 


a sound velocity 

6 airfoil chord 

pressure-drag coefficient . “ 
+1)'4Cp 


"p reduced drag coefficient (A) 


C, pressure coefficient 
or 1y^C, 
er" 


~ 
WW 


reduced pressure coefficient ( 


Gledstone-Dale constant 

model span 

Mach number 

index of refraction 

pressure 

dynamic pressure 

airfoil thickness ratio 

horizontal component perturbation velocity (per- 
turbation from a*) 

vertical component perturbation velocity 


Eaa 
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NR 
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profile 

reduced vertical distance (io + 1)t/e]*y) 

angle of attack * 

ratio of specific heats (1.4 for air) 

semiwedge angle 

wave length of monochromatic light used on inter- 
ferometer 


reduced Mach number G FIPA ) 
density 
ubscripts and superscripts: 
conditions in free stream 
reservoir conditions 
reservoir conditions behind a shock wave 
conditions at sonic velocity 
Symbois used without subscripts indicate local conditions. 
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APPARATUS AND METHODS 
WIND TUNNEL 


The measurements were made in the GALCIT 4- by 
10-inch transonic wind tunnel. For a description of the 
tunnel end the flexible nozzle employed see reference 26. 
The tunnel cen be run at both subsonic and low-supersonic 
velocities with continuous Mech number variation through 
use of the flexible nozzle and a variable second-throat nozzle 
downstream of the test section. 


MODELS 


The models used were half airfoils followed by straight 
sections. Four of the models were wedges (semiangles 
4.53°, 7.56°, 10.00°, and 26.57?) followed by straight sections 
and the fifth was half of a biconvex circular-arc airfoil (8.80 
percent thick) followed by a straight section (see fig. 1). 
The distance from the leading edge to the point where the 
straight section began was of the order of 4 inch for all five 
models. The models were made of tool steel and were very 
carefully machined and lapped so as to give exact cylindrical 


Cartesian coordinates, origin at leading edge of. 





FIGURE 1.—Geometry of sections tested, te ata lo tacha. 


surfaces. Two pressure orifices on opposite sides of the air- 
foil were placed exactly the same distance from the leading 
edge to aid in setting the model to zero angle of attack by 
balancing these pressures on an alcohol U-tube. Because 
of the very short chord lengths vernier-protractor measure- 
ments of the opening angles of the leading edge were of 
doubtful accuracy, so the angles were measured by letting 
the leading edge split a beam of parallel light and measuring 
the position of the reflected spots on: a wall behind the 


model. In this manner the angles could be measured to 


+0.03°. 
INTERFEROMETER 


The interferometer used i in this investigation is described. 
in references 27 and 28. One of the main features of this. 


interferometer is that both light beams are passed through 
the teat section, one over the model and the other ahead of 
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the model in the uniform flow field, that is, where the velocity 
is nearly the free-stream velocity. The advantages of this 
are: (1) The fringe shifts are in relation to the free-stream 
density and (2) the effects of the side-wall boundary layers 
are approximately canceled out since both beams traverse 
nearly the same boundary laver at each side window. This 
leads to improved accuracy when the interferograms are 
evaluated on the basis of the absolute value of the fringe 
shift from no-flow conditions. For these tests finite-fringe 


interferograms were used and another method of evaluation 


was devised which ts much simpler and more accurate than 
the above-mentioned technique. Infinite-fringe interfero- 
grams, while they give the constant-density contours im- 
mediately, are less accurate than the superimposed finite- 
fringe interferograms because any optical inaccuracies in the 
system cause the contour fringes to. be distorted. These 
inaccuracies are calibrated out in the superimposed finite- 
fringe interferograms. Also there are times when one does 
not know whether the density inerement between contours 
of an infinite-fringe interferogram is positive or negative; 
this trouble does not arise with the finite-fringe interfero- 
grams. A typical finite-fringe interferogram is shown in 
figure 2. 
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FIGURE 2.—Typical finite-fringe interferogram.  8.&-percent circular-arc section at 34,, 1.200. 
METHOD OF EVALUATION OF INTERFEROGRAMS 

The method of evaluation used here depends on two 
techniques: (8) Photographie superposition of disturbed 
and undisturbed interferograms and (b) fringe identification 
by a pressure measured on the model. 


REPORT 1094—-NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Direct photographie superposition of a ''no-flow" finite- 
fringe interferogram on a “with-flow” finite-fringe interfero- 
gram gives rise to dashed shadowy lines (the dashes being 
where the dark fringes of one picture cross the light fringes 
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FIGURE 3.—Typleal superimposed finite-fringe interferogram. 10% semiangle wedge at 
M y 142178. 


of the other); see figure 3 for an example of this type of 
picture. These shadowy lines can easily be shown to be lines 
of constant density for two-dimensional flow and are the 
same contours as would be obtained on an infinite-fringe 
interferogram made with perfect optical surfnees. The 
increment in density between these shadowy lines is a 
constant dependent only on the span of the model and the 
wave length of the monochromatic light being used. This is 
easily shown since the difference in optical path lengths of 
the light rays between two adjacent constant-donsity con- 
tours must be 1 wave length of the light being used. For 
two-dimensional flow the difference in optical path length 
will simply be (An, where / is the span of the model and An is 
the difference ih index of refraction between the two light 
paths. Thus 

IAn—y | (1) 


But the relation between index of refraction and density in a 
gas is given by 
n—l=kp (2) 


where k is the Gladstone-Dale constant (a function of the 
light frequency and type of gas). Therefore 
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Ap MEI (3) 


where Ap is the difference in density between two adjacent 
constant-density contours. For these experiments 


A=5461 A (mercury green line) 
k 0.1162 cu ft/slug 
[23.50 in. 


Apfp.=0.0250 per fringe shift 


where p,— 0.00211 slug per cubic foot was the usual tunnel 
stagnation density. 

The adventage of photographic superposition is not only 
in time saved but also in increased accuracy of evaluation. 
Any alight changes in fringe spacing or fringe orientation 
with respect to the no-flow interferogram which occur before 
the with-flow interferogram is taken can be almost exactly 
canceled out by causing the two superimposed interfero- 
grams to coincide exactly in & region where it is known that 
the flow was uniform, since in such regions there should be 
no isopycnic contours. This is particularly easy to do for 
supersonic flow if & portion of the flow field ahead of the nose 
shock wave is included in the interferogram. For subsonic 
flow care must be taken to include enough of the flow field 
ghead of the model in the interferogram to have some of the 
nearly undisturbed flow field for comparison; this was quite 
simple to do for the small, thin models used in these testa. 

The actual superposition technique used here was first to 
make a print (3% times enlarged) of the with-flow interfero- 
gram. This print was then placed under the enlarger and 
the no-flow interferogram negative was put into the enlarger. 
By changing the enlargement scale and moving the with-flow 
interferogram under the enlarger the fringes were made to 
coincide exactly in the regions of uniform flow. The con- 
stant-density contours could then be drawn in on the print. 
Alternatively, the first print could be made on transparent 
paper (Ansco Reprolith Ortho was used) end when the super- 
position was accomplished a piece of photosensitive paper 
was slipped under the transparent print and a print of the 
two interferograms was obtained. This was the technique 
used for figure 3. 

In order to identify the density values with the fringes a 
pressure tap was placed on each model approximately half- 
way from the leading edge to the shoulder (a region where 
the pressure gradient was expected to be large). From the 
pressure reading the density at the pressure tap wes calcu- 
lated using reservoir fluid properties (taking into account 
entropy changes through shock waves). The pressure tap 
will always lie between two fringe contours or on a contour, 
so that. by knowing the density increment between fringe 
contours, the values of the density on the adjacent contours 
ràn be obtained by interpolation. The whole interferogram 
is determined once the densitv is known on one contour (ex- 
cept for the shifts through shock waves). 

INTERFEROMETER SENSITIVITY 

It is interesting to note that the interferometric method 
has its greatest sensitivity in the transonic range. As 
pointed out previously, the density increment between two 
adjacent contour fringes is a constant 
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Now in any part of the flow field where the stagnation 

density is constant along a streamline, 
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Hence the increment in Mach number between adjacent | 


contour fringes is given approximately by 
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This function has a minimum at Ma, which is 


M=0.914 for air (y=1.4). A graph of this function is 
shown in figure 4. Note that the Mach number increment 
per fringe for these tests was always closely equal to 0.05. 
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FUTURE £.—Inerement tn Mach number per fringe contour against local Maoh number for 
alr (y=1,4). For these tests A/kip.™0.0254-0.001. 


Similarly, the expression for the increment in pressure 
coefficient between adjacent contour fringes is approxi- 
mately 


(7) 


For values of Af close to LI, this expression has a minimum 
at M.— 5 which is 1.832 for air. | 





SIDE-WALL BOUNDARY-LAYER EFFECT ON APPROXIMATING 
TWO-DIMENSIONAL FLOW 
A close approximation to two-dimensional flow over the 
whole span of the model was required since the interferometer 
integrates the value of the density from wall to wall. Ina 
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nonviscous fluid letting the model extend from wall to wall 
would theoretically give two-dimensional flow over the whole 
span. If the model did not span the whole tunnel, the flow 
would correspond to that past & model of infinite span with 
periodic gaps in it where the gaps were equal to twice the 
distance from the edge of the model to the wall. The effect 
of the side-wall boundary layers, for & model that does not 
span the tunnel, is roughly to decrease the size of this gap. 
Approximately, the gap size would be decreased by twice 
the displacement thickness of the wall boundary layer. By 
making the gap between the edge of the model and the wall 
approximately equal to the wall-boundary-layer displace- 
ment thickness, one might hope to approximate closely two- 
dimensional flow over the span. "This phenomenon is, of 
course, very much more complicated than this, particularly 
in the supersonic case where the shock waves interact with 
the wall boundary layer. However, by taking circular- 
cylinder and wedge models and varying the gap size in in- 
crements of Ne inch, it was found that the detached bow wave 
became closely two-dimensional when the gap size was X inch 
(i. e., there was no blur ahead of or behind the shock pictures) 
which is almost exactly the boundary-layer displacement 
thickness when measured without a model in the test section, 
When the gap was Ne inch the shock was blurred ahead of 
the main shock and when the gap was Ne inch it was blurred 
behind the main shock. These tests were further sub- 
stantiated by some schlieren pictures, which Mr. Walter G. 
Vincenti of the NACA Ames Aeronautical Laboratory 
kindly made available, showing a view looking down on a 
wedge model so that the leading edge of the detached shock 
appeared as a line; by varying the model span a discrete 
value of the span was found where this line was almost exactly 
parallel to the leading edge of the model, while for just 
slight variations from this gap size the shock was curved 
forward or backward. Figure 2 shows a finite-fringe inter- 
ferogram of the circular-arc section with a detached shock 
where the definition of the shock wave was unusually sharp. 
This is strong, but, of course, not conclusive, evidence that 
the flow was closely two-dimensional over most of the span. 
Further evidence that the flow differed from two-dimensional 
flow only slightly is given in the next section. 
SIDE-W ALL BOUNDARY-LAYER EFFECT ON INTERFEROGRAM 
EVALUATIONS 

À result of the method of interferogram evaluation de- 
seribed above is that the effect of the side-wall boundary Iayer 
is approximately canceled out, since the over-all fringe shift 
from no-flow conditions is unimportant, only the relative 
fringe shifts from a point of known density being used. This 


is strictly true only if the integrated side-wall boundary-layer: 


density, defined by 
T 
J pdy | (8) 


Q 


where y is the direction perpendicular to the tunnel wall and 
y=0 is the wall, is the same over the entire field of view of 


the interferometer. Obviously, this can never be exactly ` 


true since the pressure field caused by the model, the 
boundary-layer growth, and the shock-wave and boundary- 
laver interaction all tend to change this value. An indica- 
tion that all these effects might be small was obtained from 
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the model tests where pressures were measured at two 
points on the model in the center of the span, where the 
flow is closely two-dimensional; the density increment be- 
tween these two points on the model was compared with 
the density increment given by the interferogram. The 
standard deviation from zero of the difference between these 
two increments over the whole range of test Mach numbers 
was about 1 percent of the stagnation density. Also, the 
values of pressure-drag coefficient obtained interferomot- 
rically for the attached-shock-wave cases checked the 
oblique-shock theory very closely, and it is well-known that 
the oblique-shock theory checks experiment quite well. 
DETERMINATION OF FREE-STREAM MACH NUMBER 


An interesting result of the method of cvaluation just 
described is that the free-stream Mach number in subsonic 
flow can be determined from the interferogram and the 
measured pressure on the model, provided a large enough 
field of view ahead of the model is obtained in the inter- 
ferogram. This can be done by noticing that a certain 
number of compression contours appear around the leading 
edge and then expansion contours follow these toward the 
back part of the airfoil; the center fringe corresponding 
to free-stream density can then be traced out into the flow 
field (see, e. g., figs. 9(a) to 9(d) of the 10° wedge in subsonic 
flow). The exact value of the density can be determined on 
this fringe as described previously and, henec, knowing 
the stagnation density in the settling chamber, the effective 
free-stream Mach number can be determined from thie 
isentropic-flow relations. It is believed that this effective 
Mach number is & good approximation to the free-flight 
free-stream Mach number and would give the same flow as 
that measured in the wind tunnel for the very small models 
used in these tests. 

This method is more accurate at high-subsonic speeds 
than at low speeds since more contour lines are obtained on 
the airfoil at the higher speeds (see above discussion). The 
estimated accuracy in determining free-stream Mach num- 
ber in this way was +0.01 for the range of subsonic Mach 
numbers tested. 

The free stream Mach numbers for the supersonic tests 
were obtained by calibrating the flexible-nozzle jack settings 
against Mach number with a static-pressure probe in the 
center of the tunnel. The probe was traversed upstream 
and downstream in the region where the models were to be 
tested and an average Mach number was obtained there. 
The standard deviations from this average value were of the “ 
order of +0.005 in Mach number for the range of supersonic 
Mach numbers tested. 

| WIND-TUNNEL CHOKING 


In all the subsonic testing the embedded supersonic zone 
was not allowed to touch the upper or lower walls. In one 
or two of the low-supersonic tests there was a question 
whether the embedded subsonic zone touched the ceiling or 
not. In case it did, it is well-known that in such cases the 
detached shock changes its curvature near the ceiling so as 
to come in nearly normal to the walls. Sinee the models 
were so very small (hs in. thick compared with the 10-in. 
height of the tunnel), it is believed that the effect of this on 
the pressure distribution was negligible. 
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REYNOLDS NUMBER 

The value oí the Reynolds number for all of these tests 
was approximately 60,000 based on the chord of the model. 
The boundary layer on the models was laminar and no 
effort was made to trip the boundary layer to make it 
turbulent. The compression region in the shocks shown m 
the high-subsonic-flow interferograms is believed to be 
associated with the laminar boundary layer, as mentioned 
previously. 

THEORETICAL WORK ON TRANSONIC FLOW 
RELAXATION CALCULATIONS 

In 1946 Maccoll presented a paper at the Sixth Interna- 
tional Congress for Applied Mechanics in which he described 
a relaxation calculation of the compressible fow past a 20? 
semiangle wedge followed by a straight section at Mach 
numbers of 0.7 and 1.5. The flow field in both cases con- 
tained both subsonic and supersonic velocities. His main 
assumptions were: (1) Sonic velócity occurs at the shoulder 
and (2) the streamlines of the flow are perpendicular to the 
sonic line (i. e., the line where sonic velocity occurs in the 
flow). The first assumption can be shown to be correct (see 
reference 22) so that, indeed, it is not an assumption. The 
second assumption, as Maccoll realized, was only &pproxi- 
mately correct for Àf,—1.5 and certainly quite incorrect 
far away from the wedge at M, —0.7 (since the assumption 
leads to an infinite supersonic region above the wedge). In 
effect, his solution at 1f, —0.7 was “choked” in the sense 
that the back part of the body could have no influence on 
the front part. It is well-known that for bodies at high- 
subsonic speeds a finite, closed supersonic region occurs in 
the flow, so that the sonic line makes all angles possible with 
the streamlines, including 0°. The method of solution used 
was to assume positions of the shock wave and sonic line, 
calculate the residues in the relaxation net using the isen- 
tropic-flow equations (an approximation since flow behind a 
curred shock is not isentropic), and then readjust the shock- 
wave and sonic-line location, calculate again, and so forth, 
iterating until the solution closely repeated itself. Maccoll 
found that the p/p,’ distribution on the wedge surface at 
M,—1.5 was nearly identical with the p/p, distribution at 
M,-—0.7. This led him to propose that the pressure in the 
transonic region, on bodies with distinet corners, varied as 
the stagnation pressure and he presented a drag curve 
through M,,—1 for the 20° semiangle wedge calculated on 
this basis. 

Drougge in 1948, following Maccoll, calculated the flow 
past a finite cone of 45° semiangle with detached shock wave 
at M,—1.80 and Jf,—2.15, using the same assumptions 
as Maccoll (reference 18). He also made experiments on 
this cone and found the agreement with his theory rather 
good. He made several tests at lower supersonic Mach 
numbers also and found that the p/p,’ distribution on the 
cone surface did remain nearly constant except as the Mach 
number became close to the attachment Mach number. 

Drebinger in 1950 showed how to calculate, by relaxation 
techniques, the flow past finite cones and wedges with de- 
tached shocks, eliminating the isentropic-flow assumption 
and the assumption on the streamlines being perpendicular 
to the sonic line (reference 19): He calculated a specific 


example—a 26.6° semiangle wedge at Af,—1.440—and 
checked the calculated shock-wave shape and position exper- 
imentally. His calculations showed that, even for the de- 


tached-shock case, the streamlines differed from being per- 


pendicular to the sonic lines by angles as large as 30°. His 


calculation was checked in detail experimentally in these. BR 


tests and agreement was found to be excellent. 
TRANSONIC PERTURBATION THEORY 


By assuming that the velocity component parallel to the 
free-stream direction differs only by a small quantity u from 
a*, the critical velocity, and keeping only the highest-order 
terms in the differential equation, the equations of two- 
dimensional irrotational fluid motion are reduced to 


u Ou, Op 
Ot) CC 
du Ov 9) 
a” 


It was from these equations that Yon Kármán and Guderley 
independently arrived at the transonic sfmilarity laws 
(references 9 and 4). For two-dimensional steady flow past 
sections whose shape functions are the same, these laws 
imply that 
| M*— At 

- En Dice (TEX D. a d 


where 1f is the local Mach number on the surface of thie: 


(10) 


section. The similarity in pressure and drag coefficients ^ 


is then 
GC, MIA | 
CONS ln oye ee 
FDO) ,( M.—1 
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These quantities will be called reduced local Mach number, ` 


reduced free-stream Mach number, reduced pressure coef- 
ficient, and reduced drag coefficient, respectively, using 
symbols £, £,, Cy, and Cp. 

By interchanging dependent and independent variables 
in the perturbation equations, the problem becomes linear: 


(13) 


where , 


and, by eliminating x by differentiation, the Tricomi equa- 


tion is obtained: 


The main difficulties with this hodograph (u, v) plane are: - 


(14) 


(a) The mapping of physical boundaries into the hodograph 


hal 
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plane is, in general, not known until the solution to the 
problem is known so that it is not known where to apply the 
boundary conditions in the hodograph plane and (b) the 
mapping is often multivalued, complicating the solution. 
Two interesting cases are known where these difficulties 
are avoided. They are: (a) The free jet, studied by Tschap- 
lygin in 1905, and (b) the finite wedge, studied recently by 
Guderley and Yoshihara, Vincenti and Wagoner, and Cole. 
These latter studies came to the author’s attention after the 
present experimental study of the finite wedge in transonic 
flow had begun and served to make the study more interest- 
ing since the data could then be compared with the theo- 
retical results. 


THEORETICAL STUDIES OF TRANSONIC FLOW PAST THIN WEDGE SECTIONS 


Guderley was the first to formulate the problem of the 
thin finite wedge in the hodograph; he and Yoshihara found 
an approximate solution to the problem of the flow past a 
thin double-wedge profile at. zero angle of attack at Mach 
number 1 using the transonic perturbation equations 
(reference 20).. , 

Vincenti and Wagoner considered the thin double-wedge 
profile at zero angle of attack for low-supersonic Mach 
numbers where the shock wave is detached (reference 21). 
Their solutions were effected by relaxation calculations in 
the hodograph plane. Here the bow shock wave and the 
sonic line are fixed boundaries (their positions are not 
known originally in the physical plane) and the boundary 
condition on the shock is the slope of the streamlines (or 
the lines y=Constant). This boundary condition was 
first shown by Busemann, who aptly called the configura- 
tion a “hedge hog.” 

Cole (reference 22) has recently given a simple approxi- 
mate analytical solution to the flow past a thin symmetrical 


wedge followed by a straight section at high-subsonic ` 


speeds (14, €1). His solution satisfies the Tricomi equa- 
tion and the boundary conditions on the wedge and at infinity 
but not the boundary conditions on the sonic line. Effec- 
tively, his solution gives a finite vertical sonic line from the 
shoulder which is also a limiting line. Cole has indicated 
that this solution is the singular part of the solution in the 
hodograph and as such is most likely the main part of the 
solution. It is interesting to note that the drag-curve slope 
and curvature at Af_=1 obtained from Cole’s solution 
agree exactly with the values obtained from the simple 
physical considerations of the next section. Also, the pres- 
sure distribution on the wedge at M,=1 agrees within 1 
or 2 percent with that obtained by Guderley and Yoshihara. 

Since the back half of a double-wedge profile has only a 
very weak influence on the pressure distribution on the 
front half for A. 1 (only through the “last Mach wave" 
from the shoulder point to the sonic point on the detached 
shock), it is reasonable to take the solution of the double 
wedge at M, z 1 and use the front-half solutions in con- 
nection with Cole's results for Af. <1 for the wedge fol- 
lowed by a straight section and thus have a solution for the 
latter semi-infinite body completely through the transonic 
range. By using linearized subsonic theory and the shock- 
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expansion supersonic theory, the zero-angle-of-attack flow 
is obtained for all possible values of M. 

Tsien and Baron (reference 29) have shown that the 
shock-expansion theory can be expressed in the transonie 
similarity form for thin bodies in pure supersonic flow near 
M.=1. 

Von Kármán Gate 9) has indicated also how linear- 
ized subsonic- and supersonic-flow results may be written 
in the transonic similarity form since, from the Prandtl- 
Glauert similarity, in linearized subsonic theory, 





ES pa IMs d (15) 
and, in linearized supersonic theory, | 
e . & FAf fT 
Al A d (16) 
and from the expressions for reduced pressure coefficient 
. 1/3 
and Mach number, multiplying both sides by m 
these equations may also be written as: 
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y= [y+ 1e] y 


so equations (15) and (16) may be written in transonic fornt 


GE y E =E) (17) 
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The subsonic pressure-distribution and drag-coefficient 
curves have been calculated here from Cole’s analytical 
expressions and, combined with the results of Guderley and 
Yoshihara, Vincenti and Wagoner, and Tsien and Baron, 
the curves for reduced pressure and Mach number distribu- 
tion and reduced drag coefficient ? are given in figures 5 to $ 
for the finite wedge followed by a straight section. 


1 The reduced drag coefficient given In figure 8 ls that for the half wedge and ís equal to ` 


pec | C. d (5) 
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Fi;URE 5.—Theoretical reduced local Mach number distributions on a wedge near Mach 
number 1. Data from references 20 to 22 and 29. 


FIGURE 7.—Theoretical reduced pressure-coeficlent distributions on a wedge near Mach 
number 1, Data from references 20 to 22 and 29. 
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FIGUEX 8.—Theoretical reduced drag coefficlent against SE Mach number for a wedge. 





H can be shown that Cole’s solution for large negative - 
values of £, goes over exactly into the linearized subsonic 
solution (see appendix A). The reduced-pressure-coefcient 
curve for £,— — 2.02 in figure 7 is so nearly identical for both 

«Solutions that they cannot be told apart (except that Cole’s 
solution goes to Cy at rfc=1 while the linearized solution 
goes to — œ). This is to be expected since the transonic 
perturbation equations are not restricted to transonic ` 

' flow but apply equally well to completely subsonic and 
completely supersonic flow.* The transonic equation can 
be written in the form 





FIGURE 6.—Theoretical reduced local Mach number against reduced free-stream Mach 3 This was pointed out to the author by Dr. Milton Van Dyke of the NACA Ames Aerg- 
number at several choriwise stations on a wedge. Data from references 20 to 22 and 20. nautica? Laboratory. 
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where e is the perturbation potential such that u=U +3, 


LA 


Thus it is clear that for completely subsonic or 


completely supersonic flows the term on the right is negli- 
gibly small but becomes of paramount importance in transonic 
flow. 


CHARACTERISTIC FEATURES OF TRANSONIC FLOW PAST 
WEDGE AND CIRCULAR-ARC SECTIONS 


CHARACTERISTIC FREE-STREAM MACH NUMBERS 


Critical Mach number.—The Mach number at which 
sonic velocity first appears on the wedge is M,—0 (within 
the inviscid theory) since subsonic flow cannot turn a sharp 
corner. Because of the fact that the boundary layer rounds 
off the corner, and perhaps also because of the spatial 
resolution limitations of the interferometric method, sonic 
velocity was not found there experimentally until 
approximately £.—— 0.80 for the wedges. 

The critical Mach number for a half circular-are airfoil 
followed by a straight section can be obtained approximately 
from linearized subsonic theory. This theory gives the 
surface pressure distribution as 


_ — 4l) ` 


RTE MS 1-1) oe 675] eo 


which yields 


— —1.626(¢/c), | | 
Pain ar SE (21) 


at z/c—0.783 (see appendix B). This equation can also be 
written in transonic similarity form by multiplying both 


1/3 
sides by ae (as shown in the previous section): 


_ —1.626 . 


G = - . (22) 
Prin re E 
Now, within the transonic perturbation theory, 
C, =—2(E— E.) (23) 
Hence 
Ua 2f. | (24) 
Equating G Lo Cy. one obtains the critical reduced 


Mach number 
En, —0.871 


For the thickness ratio t/c=0.088 used in these tests, this 
predicts a critical Mach number of 0.834 at 2/c=0.783. 
Experimentally, the critical Mach number was found to be 
0.825 and occurred somewhere between z/e=0.75 and 0.95 
(the pressure distribution was very flat in this range). It is 
interesting to note that the experimental A. . was higher 
for the wedges than for the circular-arc profile of the same 
thickness ratio. This was probably due to a combination 
of three effects: (1) The boundary layer for the same 


. zones near the shoulder. 
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Reynolds numbers used here was fairly thick in comparison 
with the dimensions of the model and thus it “rounded off" 
the shoulder more than would be the case at higher Reynolds 
numbers. (2) The height of the supersonic zone, even for 
&n ideal nonviscous flow past thin wedges, appears to be 
quite. small until the free-stream Mach number is quite 
close to 1. This is apparent from Cole's theory and also 
from the argument in reference 25 that the height of shocks 
in the supersonie zone must be of the form 


Inf 0» 


(3) The spatial resolution of the interferometric method may 
not have been sufficient to detect very small supersonic 
There is also a large refraction 
error near the shoulder due to the high density gradients 
which tends to obscure details of the flow there. 
Shock-attachment Mach number.—The shock-attachment 
Mach number depends only on the opening angle of the 
profile at the leading edge and can be predicted quite pre- 


` cisely by oblique-shock theory. If 0 is the semiopening angle, 


then it can be shown that approximately, for thin profiles, 


EP CO 


(see appendix C). If t/c is the thickness ratio of the circular- 
arc section, 0:52(t/c). Hence for the circular-arc profile 
3 | 
Eaa 337 (26) 


Mach number at which sonic velocity appears behind an 
oblique shock.—The Mach number at which sonic velocity 
appears behind an oblique shock AM. is just slightly higher 
than Ma, and again is a function only of the opening angle. 
These values can also be found quite precisely from oblique- 
shock theory and approximately in aid form can be 
given as. e - 


M.~—1 - 
eet i, c 69 

for the wedge (see appendix C) and 
. En =2 l CLIP, (28) 


for the circular-are section. 
CHARACTERISTIC YALUES OF LOCAL MACH NUMBER 


Mach number at leading edge.— The Mach number at the 
leading. edge is zero (a stagnation point) for all free-stream 
Mach numbers less than the attachment Mach number. 

Mach number at shoulder of wedge.- The Mach number 
at the shoulder of the wedge just before the turn is always 1. 
This is easily seen in the case of flow with detached shock 
since.the only characteristic distance of the finite wedge is 
the distance from the leading edge to the shoulder which 
must determine the shock-detachment djstance, and, if the 
sonic point occurred ahead of the shoulder, the shoulder 
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This follows from the 


‘ould not influence the shock position. Subsonic flow cannot 
“urn a sharp corner so the flow must therefore reach Mech 
1umber i right at the corner. In the case of subsonic free- 
stream flow the argument is not so simple (see referenee 22). 
At the shoulder the flow around the corner is locally a 
centered Prandtl-Meyer fan starting from M—1. The Mach 
number just behind the corner is thus determined only be 
the wedge angle and is independent of the free-stream Mach 
number. Behind (his point the flow will recompress to the 
free-stream Mach number through a shock or series of shocks, 
for free-stream Mach numbers less than the attachment 
Mach number. The expression for Mach number sv 
behind an expansion from Af=1 through an angle 0 is 


21 T Af pu — 1—tan”? A < Alpe —1 (29) 


Expanding the right-hand side in powers of +/4fp;?—1 
(assumed small), the first nonzero term yields 


2 Dien 18 


ĝ = 





which is in transonic similarity form so 
Met äi 
toga" (2) BS 


PRESSURE DISTRIBUTION ON BODIES MOVING THROUGH 
AN INFINITE FLUID AT SPEEDS NEAR MACH NUMBER 1 


STATIONARY VALUE OF LOCAL MACH NUMBER AT FREE-STREAM MACH 
NUMBER 1 

During the course of these investigations it was found 
that for the wedge and circular-are sections the local Mach 
number distributions on these sections at very high subsonic 
speeds (above Af... but below choking Mach number) 
and at very low supersonic speeds (where the detached 
shock wave was a chord length or so ahead of the section) 
were nearly identical. In trying to understand why this 
should be so, the following explanation was derived: (1) At 
low-supersonic speeds the bow shock wave is detached a 
great distance &head of the profile and & subsonic flow region 
is embedded in the flow field between the shock and the 
sonic line. The part of the shock directly ahead of the 
profile is nearly normal over quite a distance (of course, the 
slope of the shock asymptotically tends to the slope of the 
Mach wave of the free-stream flow at large distances lateral 
to the flow direction). Nagamatsu (reference 30) has pre- 
viously indicated this and points out that the flow past the 
profile should be closely approximated by assuming the 
profile is in a high-speed subsonic flow where the velocity 
distribution at infinity is slightly nonuniform, the minimum 
velocity being directly ahead of the profile and equal to the 
velocity behind the normal shock and then increasing in 
both lateral directions. (2) Now the normal shock near 
Mach. number 1 is nearly symmetrical in the sense that the 
Mach number behind the shock is just as much below 1 as 


185 
the Mach number ahead is above 1. 
normal-shock relation: 


i—Mg-— Al 


— p (32) 
L 1 iú c Mr —1) 
where M, is the Mach number ahead of the shock and Af, 
is the Mach number behind the shock; so near J4=1, 
1 —Af,’? gs 11,*— 1 (33) 
1— lhe 1f, — 1 (34) 


Therefore if M, =1-+e, where e is small, the flow past the 
profile is nearly the same as the flow past the profile at 


or 


. M, —1—e since the Mach number behind the centrel part 
It. 


of the detached shock wave is almost exactly 1—e. 
follows therefore that the local Mach number distribution 
on the profile surface must have a stationery value at 


M. —1 and furthermore vary only slowly in the neighbor- - 


hood of M, —1. Mathematically this means 


em LN 0 


(35) 


It should be noticed that this argument is based on two — 
assumptions: * (1) The detached bow wave moves very far l 


ahead of the profile as the flight Mach number decreases 
toward 1. (2) The radius of curvature of the detached bow 
wave at points directly ahead of the profile becomes ex- 
tremely large as the flight Mach number decreases toward 1. 

Examining these assumptions, it would seem that the 
same reasoning should apply to any finite three-dimensional 
body id an infinite fluid traveling at speeds near Mach 


number 1, except that now two radii of curvature at points ' 


on the detached bow wave ahead of the body must be 


assumed to become large as the flight Mach number de- `` 


creases toward 1. The detached bow wave is so far away 
from the body at speeds just slightly above Mach number 1 
that the body appeers as only a very small object in relation 
to the radii of curvature of the bow wave and, hence, it 
would appear as though the shape and attitude of the body 
could have no appreciable effect in changing the argument 
presented above. 

The reasoning should also apply to an infinite yawed 
cylinder (whose cross section may be finite or, if the angle 


of attack is 0, may extend infinitely far downstream) pro- 


vided that the Mach number considered is the component 
of the Mach number normal to the generators of the cylinder. 
These arguments are for steady-flight speeds. Large 


accelerations through sonic flight speed could conceivably ` 


modify the phenomenon. Thus it is difficult to judge 
whether or not the available flight-test data confirm the 


concept since nearly all such data come from missile tests — ` 


& It Ig believed that these are not actually aesumptions but are capable of demonstration ` 


if one aseumes a smooth variation of drag through Mach number L 


wf 
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that involved large accelerations (or decelerations) through 
sonic flight speeds. The transonic-bump tests of Weaver 
on sweptback wings (reference 1) would seem to support 
these conclusions since they show drag-coefiicient maximums 
very near Mach number 1, & necessary consequence of the 
concept for finite three-dimensional bodies and finite, 
unswept, two-dimensional bodies as will now be shown. 


SLOPE OF-PRESSURE- AND DRAG-COEFFICIENT CURVES AT Ma =1 


Equation (35) enables one to calculate the slope of the 
pressure- and drag-coefficient curves at Mach number 1 
as follows: 
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Now for a two-dimensional body the pressure-drag coefficient 
(based on the chord) is given by the contour integral ' 


— $ Ct Ade 
where 


£ unit vector in stream direction 

f unit vector normal to profile pointing outward 

ds element of length along profile contour 

so if the angle of attack is constant and Af, is changing 


(39) 
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But $i de=0 for a closed contour, so 
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For the front part of a profile (defined as that part ahead 
of the maximum thickness) the usual definition of a drag 
coefficient is 


= 
A ee l 


Gan an C,t-nde (42) 


"b 
where | means the counterclockwise line integral from the 


point of maximum thickness on the upper surface to the 
point of maximum thickness on the lower surface; thus 
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where 
i maximum thickness of profile 
a angle of attack of profile 
Similarly the drag coefficient for the rear part is 
Co, — | Coto (44) 
b 
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For the tests on wedge and circular-are sections followed 
by straight sections the concept of drag cocfficient of the 
front part of the section will often be used. 

For bodies of revolution (which include spheres, eone- 
cylinders, etc.) the pressure-drag coefficient (based on maxi- 
mum cross-sectional] area) at zero angle of attack is 


d I 2 
= Y afl 
Com] i (5) 
where 


R maximum radius of body 

L length of body . 
z distance from nose along axis 
Therefore 
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as before in the two-dimensional case. However, for front 


and back drag coefficients 
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and these differ from two-dimensional values obtained above 
in equations (43) and (45) by not involving the fineness 
ratio of the body (this is of course due to the different 
reference areas for drag coefficients). 

For the general finite three-dimensional body the pressure- 
drag coefficient is given by 
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ak ff C(t-Alde (49) 


where A is some reference area of the body and S is the 
surface of the body. It follows as it did previously that 
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SLOPE OF DRAG-COEFFICIENT CURVE AT Meel IN TRANSONIC SIMILARITY 
PARAMETEES FOR TWO-DIMENSIONAL FLOWS 
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Within the transonic approximation 
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Now 
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so 
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Similarly it is easy to show that 
dC, 
dE, to a (55) 
and 
dCo, 
alar (56) 
OTHER DATA SHOWING SLOW VARIATION OF LOCAL MACH NUMBER NEAR 
Mol 


As mentioned previously, Maccoll in 1946 had already 
proposed the slow variation of local Mach number near 
M,=1 on “bodies having distinct corners." It appears 
that this latter restriction is not necessary. Maccoil's pro- 
posal was based on rather slim evidence and it is believed 
that here, on the basis of the argument presented concerning 
the normal shock, the principle is explained more convinc- 
ingly. Also the experimental evidence given here and by 
Drougge (reference 18), Bleakney and Griffith (personal 
communication), Weaver {reference 1), and by some NACA 
reports tends to bear out the conclusions of slow variation of 
local Mach number on bodies near 4f_=1. 

This fact is sometimes slightly obscured in the NACA 
reports because pressure coeffictent was plotted instead of 
pip. or local Mach number. However, constant Mach 
number lines weré sometimes drawn in these plots and there 
the evidence shows up strongly (see, e. g., reference 31, figs. 
* to 1i, pp. 36 and 37). The relative constancy of local 


212483—54———1* 
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Mach number distribution near 1,--1 for airfoils at an 
angle of attack is also shown clearly in figures 8, 9, and 10 
of reference 32. 


ON COMPARING THEORY AND EXPERIMENT 


In references 21 and 25 discussions were presented on the 
philosophy of comparing experiments with approximate 
theories, and these discussions will not be repeated here, 
except to mention that in some of the theoretical curves 
presented here the values have been shown with & certain 
spread which ue from using & pressure coefficient equal 


to ee or a9 = TF Ze 





(the former value is the one that 


fits into transonie similarity theory; the latter value is the 
one more commonly used in perturbation analysis). 


In connection with the idea presented in reference 25 of 
extrapolating experimental data to zero thickness in order 
to compare with results from transonic perturbation analyses, 
it is interesting to note that the characteristic Mach numbers 


mentioned in the section ‘Characteristic Features of Tren- 


sonic Flow past Wedge and Circular-Arc Sections” can be 


presented in powers of the thickness of the wedge (or equiv- 


alently in powers of the wedge angle), the first term of 
which gives the transonic similarity expression; two of these 
values are | 


At 
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(See appendixes D and E.) In transonic perturbation 
theory the terms in 0 on the right-hand side are neglected. 
This can lead to fairly large errors for even moderately 
large values of € since the approach to £=0 is nonlinear and 
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———»o as fl (59) 
Judging from this one might expect that quantitative 
agreement of transonic perturbation analyses with experiment 
would not be so good. However, in comparing two similar 
shapes with only slightly different thickness ratios by 
transonic similarity considerations one would expect fairly 
good agreement. 


EXPERIMENTAL RESULTS 
FLOW FIELD NEAR 10° WEDGE 
Figure 9 shows interferograms of the flow past the 10° 


semiangle wedge for 14 Mach numbers from 0.700 to 0.892 
end 1.207 to 1.465 (the interferograms for the 4%” and 


7M? wedges were very similar and hence are not shown here). . 
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(b) Me = 0.704; Poo [pe 90,743, 


FIGURE 0.—Interferograms of flow past 10° semiangle wedge for various Mach numbers. 
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Figure 9.—Continued, 
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l (d) M 0.892; p /p.="0.891 
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FIGURE 9.—Contínued. 
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FIGURE 9.—Conclnded. 


Notice that the lines of constant density in the subsonic-flow ` 


interferograms are roughly elliptical in shape as predicted 
hy the theory (see appendixes A and B). A supersonic 
flow region was first detected between Jfa —0.700 and 0.794 
(the sonic line is shown as a dashed line in the figures) and 
a shock emanating from the corner appears in the super- 
sonic zone at M,=0.794. As the Mach number was 
increased, this zone grew larger and a shock appeared at 
the rear of it, while the shock emanating from the corner 
weakened and disappeared. This rearward shock was of 
the typical X type associated with a laminar boundary 
laver, and the interferograms clearly indicated the separation 
of the boundary layer ahead of this shock. The similarity 
between the flow field at Af,=—0.892 and at Af.-1.207 
(figs. 9 (d) and 9 (e)) is striking; the base of the rearward 
shock has moved quite far back on the wake of the blunt 
trailing edge at A£,—1.207 but in the vicinity of and ahead 
of the sonic line the two fields are nearly identical except 
for the detached shock wave which appears about 1% chord 
lengths ahead of the wedge at Af.—1.207. As the Mach 
number was increased above 1.207, the detached shock 
moved in closer to the leading edge and finally “attached” 
at a Mach number quite close to the theoretical attachment 
Mach number of 1/,—1.418. Notice that the process of 
attachment is very continuous. The effect of the boundary 
layer is quite noticeable in the last few interferograms: 
This can be roughly accounted for by considering the 


boundary layer to change the shape of the body by its 
displacement thickness and then considering & nonviscous 
flow past this revised shape. On the wedge the boundary 
layer will not grow so rapidly as on a flat plate because of the 
favorable pressure gradient and, in fact, the effect of the 
strong expansion around the corner is known to cause an 
almost complete collapse of the boundary layer there. As 
the bow shock wave gets close to attachment, the velocities 
in the subsonie region behind it are getting very close to 
sonic velocity and hence the flow in this region is very 
sensitive to any slight curvature of the ''revised shape" 
of the wedge. "This accounts for the shift of the base of the 


sonic line forward to the leading edge as the shock ap- - 


proaches attachment. The nonviscous theory would in- 
dicate that the sonic line would always begin at the corner 
and, at a Mach number just slightly above the shock- 
attachment Mach number, the whole subsonic region would 


become sonic; then, with increasing Mach number, the flow ` 


behind the shoek would be completely supersonic. As 
observed, the boundary-layer effect is to make the wedge 


have a curved surface and the sonic line actually moves - 


slowly from the corner to the nose. Even with an attached ` 


shock wave at 1f, =1.465 the flow behind the shock is not. 


quite uniform (as nonviscous theory would indicate it 
should be} because of the effective curved surface caused 
by the boundary layer. 
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LOCAL MACH NUMBER DISTRIBUTIONS ON THREE THIN WEDGES 


Figure 10 shows the variation of local Mach number 
distribution on the surfaces of the 44°, 73°, and 10° semtangle 
wedges with free-stream Mach number. This should be 
compared with figure 5 which shows the corresponding 
theoretical curves in terms of the transonie similarity param- 
eters. The general behavior of the theoretical and ex- 
perimental curves is quite definitely in good agreement. 
Particularly noteworthy is the slow variation of the local 
Mach number distribution near free-stream Mach number 1. 


PRESSURE-COEFFICIENT DISTRIBUTIONS ON THREE THIN WEDGES 


The slow variation of the Mach number distribution in 
the range near Mf,=1 is obscured when the results are 
plotted in terms of pressure coefficient, since the pressure 
coefficient changes a great deal if local Mach number is con- 
stant while the free-stream Mach number changes. A better 
parameter for presenting transonic pressure distributions 
would be p/p. (p/p; in case of a detached shock). Typical 
C, distributions are shown in figure 11 for the 7X? wedge 
(the results for the 10° and 4X? wedges were very similar 
and hence they are not presented). The pomts shown were 
where the isopycnics intersected the body in the interfero- 
grams. 
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FIGUBL 10.—-Continued. 
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FIGURE 11.—Pressure distributions for 734° semtangle wedge. 


Since for a wedge the drag coefficient is proportional to the 
average C,, the drag rise is evident in the subsonic distribu- 
tions as the point where C,—0 moves rearward with increasing 
free-stream Mach number. Linearized subsonic theory 
(which predicts C5—0) locates the C,—0 point at x/c—50 
percent. Figure 7 shows theoretical reduced C, distributions 
at various reduced free-stream Mach numbers. Again the 
qualitative agreement of these curves with experiment is 
evident. 


SHOCK-DETACHMENT DISTANCE FOE THREE THIN WEDGES 


Figure 12 shows the shock-detachment distance against 
reduced free-stream Mach number for the three thin wedges 
and includes the theoretical values from reference 21. Here 
Vincenti and Wagoner's values for £, have been multiplied 


hy Sg in order to make the trensonic perturbation 
œ Alen 0* . j 

value of detachment reduced Mach number agree with the 
value from oblique-shock theory for the 7X? wedge: The 
reason for this was discussed in the section “On Comparing 
Theory and Experiment," namely, the difficulty of comparing 


3 In terms of Mach number, for the 7145 wedge the shock theory predicts attachment at 
Ma 191.38 (Eo , 1.68), white the transontc perturbation theory predicts Ma ,=1.25 (ei 
= 1.19). 


ment. 
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FIGURE 11,.—Concluded, 


transonie perturbation theory quantitatively with experi- 
Notice how rapidly the shock wave moves away 
from the wedge as the Mach number is decreased toward 1. 


DRAG-COEFFICIENT VABIATION WITH MACH NUMBER FOR THREE THIN 
WEDGES 

It was shown in reference 25 that the viscous effects on 
the wedge tend to compensate each other at the leading 
edge and the shoulder so that the over-all pressure drag is 
nearly the same as if the flow were inviscid. Thus it would 
be expected that the pressure-drag coefficients obtained by 
integrating the experimental pressure distributions would 
check the imviscid transonic perturbation theory. The 
reduced drag coefficient used here was 


á 1 
6 |. e d (5) 
which is, in essence, the reduced drag coefficient of the 


upper (or lower) half wedge. This was done since the wedge 
model was regarded as the front half of a double-wedge 


(60) 


profile and hence the value given here is the part of the re- 
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model used here. Of course, this viewpoint is valid only for 
supersonic free-stream Mach numbers. 
Figure 13 shows the reduced drag coefficients for the three 


24 thin wedges plotted against reduced Mach number. It is 
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4 , FIGURE 13.— Experimental reduced drag coefficient against reduced Mach number on a wedge. 
seen that the results give nearly a universal curve, which 
they should if the transonic similarity law is true, but that 
o 25 24 there. are. systematic variations with wedge angle. l This is 
£s to be expected based on the discussion of the section ''On 
FIGURE 12,—Shock-detachment distance against reduced Mach number for a wedge. Comparing Theory and Experiment.” The vertical lines 


through the experimental points indicate estimated accuracy 
duced drag coefficient contributed by the front half of such. | of the data. This figure should be compared with figure 8, 
an airfoil (Cp, as in equation (42)), based on the chord of the | the theoretical reduced-drag-coefficient variation with re- 
double-wedge profile, which would be twice the chord of the | duced Mach number. It is obvious that the qualitative 
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FIGURE 14.—Drag coefficient against Mach number for 444°, 734°, and 10? semiangle wedges. Comparison of theory with experiment. 
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agreement of theory and experiment is good. In figure 14 
the theory and experiment are compared directly for the 
three thin wedges. Here the theoretical drag coefficients 
are shown with a vertical spread, the upper values for 
M. >1 corresponding to the use of the pressure coefficient 


— Hy — 

(A and the lower values, to the use of the pres- 
— 2 bere e a e 

sure coefficient C,= =a Hel. The situation is vice 


versa for A. CC From this figure it is evident that the 
transonic perturbation theory gives a good approximation 
to experiment. 


FLOW FIELD AT Mao=144 FOR THE 26.6” WEDGE 


Figure 15 shows the experimental and theoretical constant- 
velocity lines in the subsonic region behind the detached 
shock wave for a 26.57? semiangle wedge at Af, =1.44. 


j | 
d 





Theoretical 
(Drebinger) 
Exper:mental 








FIGURE 15.—Constant Mach number contours for 26.60” semiangle wedge at Af 21.410. 


The theoretical analysis was made from relaxation calcu- 
lations by Drebinger (reference 19) who used the flow equa- 
tions with entropy variation behind the shock taken into 
account. The experimental constant-velocity lines were 
determined from the isopycnic lines of the interferogram by 
taking into account the lateral stagnation-pressure gradient 
behind the eurved shock. The isopyenic lines near this 


745 


strong shock wave were probably slightly in error because of 
the “smearing out" of the pressure discontinuity across the 
shock in the side-wall boundary layers. It is seen that the 


agreement between theory and experiment on detachment ` 


disíance and constant-velocity contours near the wedge is 
e | RON 
Figure 16 shows the surface pressure distribution from 
reference 17 and the present experiments. 
that the agreement is good. 


(Drebinger) 
Experimental 





FiGURE 16.— Pressure distribution on a 26.69 semlangie wedge at M a — 1.410. 


FLOW FIELD NEAR TEE 8.8-PERCENT CIRCULAR-ARC SECTION 


Figure 17 shows interferograms of the flow past the 
8.8-percent circular-gre section for 14 Mach numbers from 
0.718 to 0.936 and 1.11 to 1.500. 

Supersonic velocity first occurred at 1M, —0.825 (see the 
section “Critical Mach number”) and in figures 17 (c) and 
17 (d) a nearly symmetric supersonic zone is shown at 
M,=0.848. No shock waves were apparent in this zone, 
although a sensitive schlieren apparatus might have shown 
some weak shocks there. At 1f_=0.890 the supersonic 
zone has grown rapidly and now terminates in the A-shock 
configuration. Further increase of the Mach number to 
Af,=0.935 (figs. 17 (e) and 17 (f)) shows the supersonic 
zone increasing laterally and the terminating shock moving 
rearward into the wake of the body. Figures 17 (e) and 
17 (f) also show the density distribution at M, 21.11 (the 
detached shock wave was just out of the field of view of the 
interferometer) and it is interesting to note the similarity 
between the flow field at 1f,=0.935 and 3f,=1.11. It ` 
would appear as though the shock terminating the super- 
sonic zone at Af,=0.935 had moved rearward to form the 
trailing-edge shock (which is actually in the wake here 
because of the blunt trailing edge) and the supersonic zone 
had grown laterally until the sonic line joined with the 
detached shock far away from the body at M,-—1, thus’ 
causing an embedded subsonic zone in the supersonic flow 
with further increase in Mach number. | 


Again it is seen 
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FURR 18.— Local Mach number against de for increasing free-stream Mach number, 83- 
percent circular-arc section. 
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FIGURE 17.—Concluded. 


With further increase of Mach number above Af. — 1.11, 


figures 17 (g) to 17 (n) show that the detached shock again 


approached the leading edge and the embedded subsonic 
zone decreased in size until finally the shock “attached”” 
somewhere between 44, =1.400 and M,-—1.480 (the theo- 
retical value being M, es 1.423). 


LOCAL MACH NUMBER DISTRIBUTIONS ON 8&.58- PERCENT CIRCULAR-ARC 
SECTION 


Figure 18 shows the local Mach number distributions for 
the 8.8-percent circuler-arc section as obtained from the 
experiments at various free-stream Mach numbers. Again 
it is apparent that the variation of local Mach number 
distribution near 4f,=1 is very slow and, indeed, the dis- 
tribution for A7. —1 could be interpolated from this figure 
with good accuracy. 

Figure 19 is a cross plot of the data of figure 18 except 





FIGURE 19.— Variation of reduced local Mach nnmber distributions with reduced free-stream 
Mach number. 8.8-percent circular-are section. l 
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that here the data are given in transonic similarity param- 
eters. This figure shows contours of constant reduced 
local Mach number on a plot of reduced free-stream Mach 
number against chordwise position. The dashed lines repre- 
sent subsonic local Mach numbers; the solid lines, supersonic 
local Mach numbers. Note again the slow variation of 
local Mach number distribution with free-stream Mach 
number near sonic velocity. | 
PRESSURE-COEFFICIENT DISTRIBUTIONS ON $.8-PERCENT CIRCULAR-ARC 
SECTION 

Figure 20 shows the pressure-coefficient distributions on 
the 8.8-percent circular-arc section for various free-stream 
Mach numbers. The points shown are where the isopycnics 
intersected the body in the interferogram. Again the 
presentation in this manner obscures the interesting fact 
observed in figure 18. | 
DRAG-COEFFICIENT VARIATION WITH MACH NUMBER FOR 8.3-PERCENT 

CIRCULAR-ARC SECTION 

Figure 21 shows the experimental determination of the 
drag coefficient of the front part of the 8.8-percent circular- 
arc section. This again is of the nature of a fore drag coef- 
ficient and, as shown in equation (43), it should have a 
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FIGURE 20.— Pressure distributions on an 8.8-percent cireular-arc section, 
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(b) Supersonic, 
FIGURE 20.—Concluded. 
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is how the subsonic data have been joined with the super- 
sonic data. The vertical lines through the experimental 
points again indicate estimated accuracy of the data. For 
the case of an attached shock the pressure distribution can 
be calculated using characteristics theory and the shock 
polar; however, a close approximation is obtained by con- 
sidering the flow behind the shock wave to be Prandtl- 
Meyer flow. (This yields, approximately, parabolic-shaped 
bow and trailing-edge shock waves; see reference 33.) From 
this pressure distribution the drag was calculated and is 
shown in figure 21. Taking into account the “reflected” 
characteristics from the shock wave would give more com- 
pression and increase the drag coefficient so that it would 
agree better with the experimental values at M. — 1.450 and 
1.500 shown in figure 21. 

Note that the tests were made at low enough supersonic 
speeds to get definitely below the drag-cocfficient maximum 
at M, 71.20. Ñ 
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FiGCRE 21.— Drag coefficient against Mach number for 8.8-percent circular-are section. 





OCAL MACH NUMBER DISTRIBUTIONS ON A 1-PERCENT BICONVEX 
CIRCULAR-ABC AIRFOIL 

Figure 22 shows local Mach number distributions from 
reference 8 for high-subsonic-speed flow over a 12-percent 
biconvex circuler-are airfoil (with turbulent boundary 
layer). The data for the 8.8-percent circular-are section at 
two supersonic speeds have been scaled according to the 
transonic similarity laws to the 12-percent case and are 
shown for the front half of the 12-percent airfoil in figure 22. 
The back half for these two cases has been faired in using a 
Prandtl-Meyer expansion which should be approximately 
correct (& more accurate determination could have been 
made using characteristics theory and the shock polar). 
At M, =1.58, the theory indicates that the shock is attached 


with sonic speed just behind the shock on the leading edge, 
so that the distribution can be obtained by standard methods ` 


mentioned above; again the Prandtl-Meyer expansion ap- 
proximation was used for the distribution at 4Af,=1.58 in 
figure 22. 

The behavior of the Mach number distributions is similar 
to that of the distributions shown previously, except in this 
case the movement of the shock terminating the local super- 
sonic zone is shown. Apparently little change in local Mach 
number distribution occurs between .M,--0.936 and 
M, e 1.29. 


DRAG-COEFFICIENT VARIATION WITH MACH NUMBER FOR A 12-PERCENT 
" BICONVEX CIBCULAB-ARC AIRFOIL 

The data of figure 22 were converted to pressures which 
were integrated to give the pressure-drag coefficient for the 
various free-stream Mach numbers. The results are shown 
in figure 23. In addition the drags of the front and back 
halves are shown separately. The drag-coefficient variation 
between 44, —0.96 and Af, =1.20 was based on constant local 
Mach number distribution at values interpolated between the 
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FIGURE 22.—Local Mach number against r/e for Increasing free-stream Mach number. 
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o Cp of whole airfoil 
9 Co of rear half of airfoil 


6 Ca of front half of airfoil 
—-- —— Prandtl- Meyer flow 
approximation 
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FIGURE 23.—Drag ooefficient against Mack number for 12-percent bleonvex clrcular-are airfoil, (From data of Ag, 22 on basis of constant Mach number distribution from M. «0.98 to 1.20.) 


curves for Af, =0.936 and M,,—1.29. The data were faired 
into the curves for attached shock wave calculated on the 
Prandtl-Meyer expansion basis. It is seen that the fore drag 
coefficient has a maximum after M,,=1 while the drag coeffi- 
cient of the rear part has a maximum before M,=1. The 
over-all airfoil has a maximum drag coefficient just before 
M.=1 in order for the curve to have the slightly negative 
slope at 414, —1 given by equation (41). 
CONCLUSIONS 


An experimental investigation of transonic flow past. two- 
dimensional wedge and circular-are sections was made using 
a Mach-Zehnder interferometer. The conclusions may be 
stated as follows: 

1. The transonic similarity theory of Von Kármán and 
Guderley was checked and found to be in good agreement 
with experiment for thin wedge profiles near a free-stream 
Mach number of 1. 

2. The results of theoretical calculations, using transonic 
perturbation theory, made by Guderley and Yoshihara, 
Vincenti and Wagoner, and Cole for a wedge in transonic 
flow were checked experimentally at high-subsonic and low- 
supersonic speeds for three wedges of different angles and 
were found to be in good agreement with experiment. 

3. The flow field and the surface pressure distribution for 
a 26.6° semiangle wedge at a free-stream Mach number of 
1.44 were obtained experimentally and were found to be in 
excellent agreement with the theoretical calculations of this 
flow made by Drebinger. 


number 1. 


4. The pressure distributions and drag coefficients for an 
8.8-percent circular-arc section followed by a straight section 
and for a 12-percent biconvex circular-arc airfoil were pre- 
sented completely through the transonic range. It was 
shown that some difficulty arises in comparing two-dimen- 
sional transonic perturbation theory with experiment, since 
this theory neglects thickness-ratio terms of order (t/c)!^? 
and higher; for even moderate thickness ratios this will 
cause noticeable deviations from more exact theory. 

5. It was shown from some physical arguments that the 
local Mach number distribution on bodies traveling through 
an infinite fluid has a stationary value at free-stream Mach 
This was verified experimentally for the case of 
two-dimensional flow. It was shown that this concept implies 
a drag-coefficient maximum just below free-stream Mach 
number 1 for all bodies in steady flight. This fact can be 
used to obtain the variation of local Mach number distribu- 
tion on bodies completely through the transonic range of 
velocities from wind-tunnel tests, provided small models 
are used so that tests can be carried well above critical Mach 
number and to low enough supersonic Mach numbers. so 
that the bow shock wave is detached a chord length or so. 


CALIFORNIA INSTITUTE OF TECHNOLOGY, 
PASADENA, CALIF., June 1, 1951. 


APPENDIX A 


ASYMPTOTIC REPRESENTATION OF COLE'S SOLUTION FOR LARGE NEGATIVE VALUES OF REDUCED FREE-STREAM 
MACH NUMBER e 


Cole’s solution for the high-subsonic-velocity flow past a, 
thin wedge (reference 22) is given as follows (in Cole’s 
notation): 


222, y vf, Sh sinh we SEET 1 (A2)AdA 
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where the center line of the wedge is at y—0; the leading 
edge, at r—0; the shoulder, at r=1; and 


S: 0-32 SEA 
230—247) =3| o EN 


i Ia 
z -$0- M e| - 03 1) Sel | (63) 
t.=(y T 1)0 


and the other notation is the same as that in the present 
paper. 
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Using the standard methods of partial-fraction expansion, 
one may write 
sinh A(r,— 
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Substituting these into the integrals above and making use 
of the "im 
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Making use of the asymptotic formulas 
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equations (61) and (62) can be written as 
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one can write’ equations (68) and (69) for large values of g 
and z, as 
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! Equation (67) is from reference 34; equation (65) is obtained by differentiation of equation (67) with respect to y. 


1 Figures 5, 6, and 7 were calculated from this equation for :=0, for the cases E, 20 
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T Ru mah 
1 (E : ‘2, >2-70° 
T sh 7S 081 1> 
co 
Ec ic i (75) 
SE 1/6 cos Y Ce Wi 
sh éa) ef 22 8 
————— 008 T — 
Us Uo 


Eliminate » between equations (74) and (75) for #1 52 œ 


2 PEG e 4 


6 y o +3 sinh 2 sinh [r(2—2)/00] 
la Lied" (leif ` 
2(82,/2)/* sinh [x(z — 23)/v,] 2 sinh [x(z— —21)/ Pol 
76) 


ine for 2>2,>0 simply replace z by z—1 and z—z by 
2 in equation (76). Thus the lines of constant Mach 
Ata are ellipses with centers on yet) with ratio of semi- 


axes equal to : | 
(22,2) ^J 1— R2 | (77) 


which is precisely the solution given by the linearized 
subsonic theory (see appendix B). 


Now, in the notation of the present paper, 
a—a) 2 
El, PE (C gu (s 


and since £— £, is small on the wedge and since 


C, —2(E— E.) 


one can write 


ru e) PE pus | E cb á St 
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Hence, 
rlż— z) 2r 2 C, | | 
EA cl, £(—&) — j 


so, approximately, 


(e le 
Do 


AS pat = ZONA RU 


Similarly, for large values of z and 2 it follows that 


1/6 
(2 £z 1; 2,2,— co for Z se äi small (79) 
1 


Substituting equations (78) and (79) into equation (16), one 


gets the exact linearized subsonic solution for constant- 
velocity lines (see appendix B). Therefore on the wedge 
(y=0), from equations (76), (78), and (79), one has approxi- 
mately for large values of z and z 


1 1 
NK eens > etal = : 
tato 7 JG ES, 


Solving this for Cy, 





or 


C. em. (80) 


which is precisely the linearized subsonic solution for flow 
past a wedge (see appendix B). Thus Cole's solution far 
away from M,=1 tends exactly to the linearized subsonic 
solution. 


APPENDIX B 
LINEARIZED SUBSONIC AND SUPERSONIC FLOW PAST WEDGE AND CIRCULAR-ARC SECTIONS 


LINEARIZED SUBSONIC FLOW PAST A WEDGE 

Let the wedge center line be on y=0, with the leading 
edge at z=0 and the shoulder at =c. Then the incom- 
pressible-flow problem is to find an analytic function u—iv 


such that v=0 on y=0 except for 0<1x<e where v— UH 


and «—iv=0 at infinity. Such a function is 


m Sie ` 
d | C log, —H G-I (81) 
where 2=x-+1y and U i is the free-stream velocity. Thus on 
y=0, 
_ wje ` 


—(z/e) 


Using the Prandtl-Glauert transformation, for linearized 
subsonic flow 


=— 16 % TLIA (82) 


O. C, 8 — 20 lo zie 
"TEMA SM 1589 


or in transonic similarity notation 


(88) 


x —2 lon zje 
E, T 1—(z/e) 
For the incompressible case the lines of constant pressure 
in the fluid will be where 








has 


(84) 


2fe 
(z/c)—1 


but these are circles with centers at 


- Constant 


e T | o 
reese, Yy=0 | (85) 


and radii + 


In the Prandtl-Glauert transformation the y distance is 
stretched by the factor /1—44,* as is the pressure coef- 
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cient so the lines of constant pressure (and hence density) 


are ellipses with ratio of axes equal to ¥ —M. given by 
the equation 


— — —— j 
is [2451—M 3 sinh (zC,/1— M.E STS 
i ¿707 TR A55 i 
"2 sinh (+01 M26) | —1 (86) 
i/sinh (x C,J1— M WG 


LINEABIZED SUPERSONIC FLOW PAST A WEDGE 


From the Ackeret theory the pressure coefficient in super- 
sonic flow is proportional to the slope and for the wedge 
yields simply 
20 


C = 87 
or 
e (88) 
Ata 


LINEARIZED SUBSONIC FLOW PAST A CIRCULAR-ARC SECTION 


For the circular-arc section, the slope of the surface 
varies almost Hnearly with distance from the zero-slope 
point along the axis of the profile. For the section shown 
in figure 1 then, with the center line on y=0, the leading 
edge at 20, and the zero-slope point at z—c, the incom- 
pressible-flow problem is again to find an analytic function 
u—iv such that on y=0, v=0 except for 0<z<e where 





v=2U « E) where t is the half thickness at *=e and 4(t/c) T 
iw! | 0) 2 ) (95) 
u—iv=0 at infinity. Such a function is 4M —1 C 
| or 
uip U t[(2— sje GG d 96 
twa — Ie 1) log, Sec (89) > VE. : (96) 
APPENDIX C 
TRANSONIC SHOCK POLAR 
'l'he equation of the shock polar in the hodograph plane is | Letting 
I œ u' —(y 1) zl 
or yy E _ T (100) 
Fu) a* y 2 U ü (97) | +0 as) 
lat ab oue. then has "LES 
Ar P= Ua — u^) (u' Ee) (101) 


where U is the velocity ahead of the shock and Y and ? are 
velocity components behind the shock parallel and perpen- 
dicular to the direction of U, respectively. Making the 
transonic approximation in this equation, let 


u=a*+u ' 
r=? (98) 
LU —a* Fu, 


Substituting into equation (97), neglecting higher powers of 
the perturbation velocities, one obtains i 


a Al 


Ea (u,—u)? (ui) (99) 
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On the wedge (y=0,0<r<e), then 


ote) lt] om 


so the linearized subsonic. solution is 








ETE "IDEE 
o alL Bas, ASTU OD 
or 
ia me eet | T 


The minimum €, is obtained by differentiation, and one 
finds that 


at the point where log, 1-45 zic 


solution of this transcendental equation is 


and numerically the 


x/c=0.783 (93) 


- which gives 


| —1.626(t/c) 


Cra RM 


(94) 


LINEABIZED SUPERSONIC FLOW PAST A CIRCULAR-ARC SECTION 


The result here is again simple from the Ackeret theory: 


The wedge angle for detachment of the shock will now 
be given by the maximum value of v’. This is easily seen 


to occur at Waa.’ giving Ones! =5 13 a HT, Since l 


within the transonic approximation 


p' —(y- 1)0 
on the wedge and | 
ag fe M *—1 (102) 
this implies that 
^ Dr Daama (OL LPS (103) 


or, viewed in another light, this implies that the reduced —. 
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attachment Mach number is 


1.19 (104) 


| M.j—1 3 
WEE BAC 
Similarly, the wedge angle for obtaining exactly sonic ve- 
locity behind the shock is given by the value of o where 


u’=0. This is v’ == (u,^)"*. Again using equation (102), 
Y 
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this can be written 


(y+ IDa = XO. SD (105) 


and, viewed in another - this implies that the reduced 


Mach number for which sonic velocity is obtained behind, 
the shock on & wedge is 


M..— 
E bai E Ze 1.26 


Laer TS (907 


APPENDIX D 


VARIATION OF REDUCED MACH NUMBER AT WHICH SONIC VELOCITY OCCURS BEHIND AN OBLIQUE SHOCK WITH 
FLOW-DEFLECTION ANGLE 


The oblique-shock relations can be written 


| 14754 My 





+ M, cost 8 
yM? sin? pt a jit. iM sin? 


My sin” E 1 
tan 6=2 cot Bx y Ecos ES) 


JA 107) 


(108) 


where 


M: Mach number ahead of shock 
Ma Mach number behind shock 
B shock-wave angle 

6 flow deflection angle 








where 


jy i e a td E e. 


Ha- 
Expangding the right-hand side in powers of M,?—1 (assumed 


small), one obtains 


1+107 
8(y-F 1) 


Reverting this series and letting Zeien @ and Af,—M 
one finds 


M eg 7 1 
fs o DS ek 





(+1) tan e 057 ^T , Me- p...) (110) 


ag; 


—_1+10y ` le 8+ 0 ^| D 


12(75*) 


For M;—1, eliminating 8 between these two equations yields 
| ids For y=1.4, 
tan bo EAT: i^ (109) VER 1 STR = 1.176 
“iii”, 1 +3 a. AM) af 
APPENDIX E 


VARIATION OF REDUCED MACH NUMBER BEHIND A PRANDTL-MEYER EXPANSION FROM M=1 THROUGH AN 


ANGLE 0, WITH 6 


The exact relation here 1s 


e Ei tan” 1 MI i—tan- MILII 


Expanding the right-hand side in terms of 444 py?—1 , using 


q? i 


tan”! r=x A a 
one obtains 


+ 10= Mé 19^ E 


næ] 


1)*a4(M pry? — 1)* 
Tel (113) 


——- 


3 


Bux! LT e 
SE I SE 


where 


=1 Pu 


Reverting this series, the first few terms are 


Mrt [ID dp fit 1+ ln HL Dol ow) 


(114) 


TE 1) 


Therefore, 


t _ Mrl, : 
n or D 


ey po 


For y=1.4, 











4v HUT. E 
SCT — 1.09” 
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